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PsJ ' Abstract We will investigate the influence of the inhomogeneity of the universe, es- 

^ ■ pecially that of the Lemaitre-Tolman-Bondi (LTB) model, on a gravitationally bound 

' local system such as the solar system. We concentrate on the dynamical perturbation 

\^ ' to the planetary motion and derive the leading order effect generated from the LTB 

model. It will be shown that there appear not only a well-known cosmological effect 
^ arisen from the homogeneous and isotropic model, such as the Robertson- Walker 

(RW) model, but also the additional terms due to the radial inhomogeneity of the LTB 
Vh^' model. We will also apply the obtained results to the problem of secular increase in 

the astronomical unit, reported by Krasinsky and Bmmberg (2004), and imply that 
the inhomogeneity of the universe cannot have a significant effect for explaining the 
' observed dA\J/dt = 15 ±4 [m/century] 

\^ • Keywords Celestial Mechanics • Gravitation • Cosmology • LTB Model 

CO ' Ephemerides • Astronomical Unit 

^' 
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1 Introduction 

The advancements in astronomical and astrophysical measurement techniques, par- 
ticularly those involving with the solar system, has been achieved remarkable ac- 
curacy of up to 9 to 11 digits level. These technical advancements have drastically 
improved the accuracy of planetary ephemerides such as DE ||T], EPM ||2l, VSOP Ol 
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and INPOP (|4] and that of various astronomical constants. With the increasingly im- 
proved measurement techniques, observational models are also required to be more 
accurate and rigorous; for details, refer to ["5*1 and the references therein. 

High-precision observational data also play a crucial role in experimental rela- 
tivity 16il7 |. Presently, the main parameters of parametrized post-Newtonian (PPN) 
approximation, j3 and 7 are tightly constrained to the value of general relativity, i.e., 
P — y— 1. For a more accurate verification of gravity, space tests such as LISA Ii8l . 
LATOR H, and ASTROD/ASTROD i 1 10| have been planned. 

Thus far, theoretical developments in studies of the solar system have pertained to 
slow motion, slow rotation and weak field approximation, g^v ~ rj/iv+hnv, \hfiv\ ^ 
1, where Tj^v is the static Minkowski metric and h^v is the perturbation. See ifTTI 
[T2l[T3l[T4ll 1 5 |[T6l . However, it is well known that our universe is expanding at an ac- 
celerated rate IfTTI . Therefore, it is natural to consider the situation that the metric 
tensor, instead of the Minkowskian metric, asymptotically reaches for the expanding 
spacetime or the background Minkowskian metric rji^y is replaced by the cosmo- 
logical type. Several investigations have been conducted, which combine the local 
metric, e.g., the Schwarzschild spacetime or the barycentric celestial reference sys- 
tem adopted by lAU, with the global cosmological comoving coordinates Ill8lll9ll20l 
|2l]|22l|23]|24l|25l|2i|22l|28l|23 . 

The cosmological contribution to the local system has thus far been discussed 
based on the homogeneous and isotropic cosmological model, i.e., the Robertson- 
Walker (RW) model. However, inhomogeneous cosmological models have recently 
attracted considerable attention since these models can provide a possibility to expla- 
nation for the observed accelerated cosmic expansion without introducing the concept 
of dark energy. For instance, the luminosity-distance was investigated in Il33ll34ll35l 
[36ll37l based on the local void model and in fiE\ using the Lemaitre-Tolman-Bondi 
(LTB) model II39II40II41II42L Moreover, Kasai (2007) 143J re-analyzed the observed 
Type la supemovae data and proposed a phenomenological method to describe the 
large-scale inhomogeneity of the universe. 

Therefore, it would be significant and interesting to investigate the influence of 
the inhomogeneity of the universe on the gravitationally bound local system. As far 
as we know, this issue has previously been examined by Gautreau (1984) ll44l and 
Mashhoon et al. (2007) ll45l . Gautreau studied the special case of the LTB model, 
<S{r) = 0, (see ( fT9] l) with respect to his cosmological theory of the curvature coordi- 
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nateo While, Mashhoon et al. investigated the cosmological contribution due to the 
LTB model as the tidal dynamics in the Fermi normal coordinate system. 

With the remarkable improvements in the observations, it has been found that 
there exist the unexplained phenomena in theory within the solar system; the pioneer 
anomaly |47 1, the Earth fly-by anomaly pS), the secular increase in the astronomical 
unit 1,49 i . and the anomalous perihelion precession of Saturn I.50J . Presently, the ori- 
gins of these anomalies are far from clear. Nonetheless, they may be attributable to 
some fundamental properties of gravitation (see llSTI and the references therein). 

Among such phenomena, the secular increase in the astronomical unit is of con- 
cern to us. From the analysis of radiometric data, Krasinsky and Brumberg ll52l dis- 
covered the positive secular trend in AU aj^ 

dK\] , , , 

^ = 15±4[m/cy], (1) 

see also ll52l . Recently, Pitjeva and Standish evaluated dPi\J /dt ~ 20 [m/cy] f53l. 
These estimated values are approximately 100 times the error of the present best-fit 
value of AU 12], 

EE AU = 1.495978706960 X 10^' ±0.1. (2) 

1 [m] 

This secular trend in AU was found by using following relation ll54l 



_ '^theo 
'theo — 



[s], (3) 



where ftheo is the computed value of the round-trip time of the light/signal (theo- 
retical value); c/theo- the interplanetary distance evaluated from the lunar-planetary 
ephemerides in the unit of [AU]; c, the speed of light in vacuum; and fo, the initial 
epoch of ephemerides. AU and dAV/dt are, respectively, the astronomical unit and 
its time variation, ftheo is compared with the observed lapse time fobs- 

The time dependent term in ^ cannot be correlated with any theoretical predic- 
tion; hence, several attempts have been made to explain this phenomenon, such as the 
effects of the cosmological expansion 1491145 Il55l . mass loss of the Sun II49II56L the 
time variation of gravitational constant G p^, and the influence of dark matter ^57\. 
However, none of these attempts have thus far been successful. 



' Gautreau does not start from tlie original form of the LTB model. However, Krasiiiski |46 | suggested 
that the model by Gautreau corresponds to the sub case of the LTB model, iS'{r) = 0. 

^ In this paper, cy refers to century as used in the study of Krasinsky and Brumberg (2004). 
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Krasinsky and Brumberg (2004) pointed out that the inhomogeneity or non-uniformity 
of the universe may have a possible explanation for dAU/dt; however, they did not 
provide evidence to support this hypothesis. Hence, it is important to verify this indi- 
cation. Further a clarification of the observational difference between homogeneous 
and inhomogeneous cosmological models in the local dynamics is important in the 
field of modem cosmology. 

In this paper, we will focus on the LTB solution as the inhomogeneous cosmolog- 
ical model and investigate its contribution to planetary motion. In section |2] we will 
summarize the dynamical perturbation based on the isotropic and homogeneous RW 
model and Robertson-McVittie (RM) model. Next, in section[3] we will derive the dy- 
namical perturbation attributed to the LTB model. As an application of the obtained 
results, we will consider the secular increase in the astronomical unit, reported by 
Krasinsky and Brumberg (2004), in section |4] Finally in section |5] we will conclude 
the paper. 

2 Dynamical Perturbation in the RW and RM Models 

Before discussing the dynamical perturbation in the LTB model, let us provide a brief 
overview of the planetary perturbation due to the Robertson-Walker (RW) model, and 
subsequently, due to the Robertson-McVittie (RM) model. Without loss of generality, 
we first consider the flat (k = 0) RW metric in the standard comoving form, 

ds^ = -c^dt^ + a^it)[dr^ + r^{de^ + sm^ ed^^)], (4) 

where a{t) is the scale factor The equation of motion of a test particle can be ex- 
pressed as Il58ll59l . 

d'^x' dx^ dx^ 1 dx^ dx^ dx' 

'W ~ ~ ^"'dT'di^c dt dt 'dt 

where t is the coordinate time; F^y, the Christoffel symbol; and v' the coordinate 
velocity. Restricting the equatorial motion to = n/2, the equations of motion for r 
and (j) are given by 
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where H = a/a is the Hubble parameter. If we introduce the proper or radial length 
R as. 



R = m{t), 



then (|6]l and d?) are rewritten as, 



cf-R /d(j)\ 



dt I 



d_ fj^2^ 

dt V dt 



a 
a 



0. 



(8) 

(9) 
(10) 



Hence, from the point of view of R, the motion of the test particle in RW spacetime 
is governed by 

(11) 



Next, in order to observe the effect of cosmological expansion on the Newtonian 
gravity, we adopt the RM solution B60II18L 



ds-" 



1 - 



GM 

2pTa(t) 



GM 

2c^ra(t) 



c^dt^ - 



GM 



l4 



2c^ra{t) 



a^{t)idr^ + r^dQ^), (12) 



where dQ^ = dO^ + sin 9d(j)^, G is the Newtonian gravitational constant, and M 
is the mass of the central gravitating body, i.e., the Sun. (fT2l i is expressed in the 
Newtonian or 1 st order approximation as 

2GM ' 



1 



c^ra{t) 



c^dt^ + a^{t){dr^ + r'-dn^). 



(13) 



(ITTb can be alternatively obtained from the cosmological perturbation theories IISTI 
|62l|63]|Mll65]|66ll3Ql, 



ds^ = -[I +21'{t,x)ydt^ + a\t)[l +2<P{t,x)]5ijdx'dx', 



(14) 



where f relates to the Newtonian gravitational potential, is the perturbation to the 
spatial curvature, and 5ij is the Kronecker's delta symbol. 

From (|5]l and ( fT3l l. the equations of motion can be expressed as 
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Further, using dHJ, the coordinates of (flSt and (IT6^ are transformed into the proper 
coordinates as 



From ( [TtI i and (fTsT l. we find that in the Newtonian or 1st order approximation, the 
leading term of dynamical perturbation obtained from the RM model is the same as 

fRWl fRW) 

those generated from the RW model, Fg and 
3 Dynamical Perturbation in tlie LTB Model 

It is generally difficult to construct any cosmological model containing a gravitat- 
ing body because of the non-linearity of general relativity; however, the Robertson- 
McVittie (RM) model is an exception. As shown in the previous section, it may be 
a practical working hypothesis that the equation of motion due to both gravitating 
body and the cosmological effect can be phenomenologically determined by the lin- 
ear combination of the Newtonian gravitational attraction and the cosmological effect 
evaluated in the cosmological background metric without considering the gravitating 
body. 

With the above assumption, let us obtain the cosmological perturbations attributed 

(RW) fRW) 

to the LTB model, which can be used to replace Fg ,F^ given in the previous 
section. 

The metric of LTB spacetime in the standard comoving form is given by II39II40I 




(17) 



(18) 



mm. 



ds^ = -c^dp- + 



l+2S{r) \ dr 




(19) 



Here, M denotes the functions of / and r, and 




(21) 



(20) 
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in which A is the cosmological constant, p{t,r) is the density of the cosmological 



pressureless particle^ and <o'(r) and ./#(r) are the arbitrary functions of r. S{r) 
is the generahzation of the curvature parameter k in the RW model, and ./# (r) is 
the active gravitational mass that generates the gravitational field. It may be noted 
that ^ has the dimension of physical length, namely, the source area distance or the 
luminosity distance, while r is a dimensionless coordinate value|42|. 
Using ^ and ( fT9] l. we obtain the equations of motion for r and ^ as 



dt^ 



dr d^M 



' dtdr dt 

1 dS 
\+2S~d^ 



2 



dt 



dr 

dt 
dM dr 
dr dt 



d^ 
'dt' 



d(j) 
dt 



1 

"37 



(22) 
(23) 



where we ignored the ff{c ^) and higher order terms. When we put the flat RW limit 
m~^R^ ra{t),S -> ;t = 0, ^ and (|23]l reduce to ^ and ([TO]l, respectively. 

In order to relate r to ^ explicitly, we suppose that the background LTB spacetime 
is regular at the origin r = where the central body is located, and that the test 
particle, such as a planet, moves around r = 0; hence, the cosmological redshift z in 
this area is sufficiently small, z ^ 1. Thus according to fl31 . we adopt the following 
expansion forms for S" and ^ around r = as 



S^{t,r) = ra{t) 



1 



1 



-\±-A{t)r+ff{r^) 
2 ayt) 



A{t) = 



dr^ 



«1, 



S{r) = -£r^ + , 



r) = —mr 
6 



(r3), e = 



dr^ 
d^yi 



dr^ 



«1, 



<1, 



in which the scale factor a{t) is defined as, 

, , dm 

«(0 - ^ 



r=0 



Using these equations, (l22l i and (l23l l are rewritten as 

/ja\2 



d(^ 
dt 



(LTB) 



^ In tlie case of tlie LTB model, the energy-momentum tensor is given by T^^ 
p(t,r) is the density and is the 4- velocity. 



(24) 
(25) 
(26) 

(27) 



(28) 

- p{t ,r)u^u^ , where 
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dt 



dt 



LTB) 



(29) 



where the leading-order dynamical perturbations, ,!^!t™^ and ^1^^^' are expressed 



as 



^(LTB) _ 



1 dA 
A dt 



2e 





n dA'' 




^ 2e 


^fl2 ^ 




\A dt ^ 


)] 


^ 

a'^ a 



In (|30] |. we used the standard relation in the RW model, 

2 



a 



(30) 
(31) 

(32) 



where q is the deceleration parameter The first term in (1301 is F^'^\ and the second 
to fourth terms are corrections obtained from the LTB model. It may be considered 
that the second term in ( l30l l is analogous to F^*^*' ~ ~q{d/a)^. 

In (l30l l. we must evaluate e and A{t). Since the observational cosmology indicates 
that our universe has a flat geometry, we can set e = 0. 4(f) may be in principle 
obtained from the modified luminosity -redshift relation [67], 

1 dA 



d\. = c 



iH^ dt 



(33) 



Finally, following the assumption discussed in the beginning of this section, the equa- 
tion of motion, attributed to both the gravitating body and the cosmological effect due 
to the LTB model, can be phenomenologically given by, 

d^^ ^(d^^^ 
dt^ V dt 



GM 



dt \ dt 



(34) 
(35) 



4 Secular Increase in Astronomical Unit 

In this section, as an application of ( |34] | and ( l35b . let us consider the secular increase 
in the astronomical unit 049ll52ll55l . Krasinsky and Brumberg found from their anal- 
ysis of planetary radar and martian orbiters/landers range data that the astronomical 
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unit (AU) increases with respect to meters as dAU/dt = 15 ±4 [m/cy]. This secular 
trend cannot be related to any theoretical model, and thus far, the origin of this secular 
increase is far from clear 

Krasinsky and Brumberg suggested one possibility, that the inhomogeneity of the 
universe may be an explanation for dAU/dt. We consider this possibility in terms of 
the LTB model. Because the current cosmological observations assert that the geome- 
try of our Universe is flat, we choose e = 0. In this case, the cosmological perturbation 
=^1^^™' becomes 



<y(LTB) 



flV (\dA 
a) \A dt 



(36) 



In our approximation, —q{d/ df^M or, equivalently, {d/a)M is a dominant cosmo- 
logical effect; however, its contribution is considerably negligible 111911201121112511261 
I32ll27ll29ll28l[30ll55l . From the assumption ( l24l i. A can be considered as a small cor- 
rection to the scale factor a{t); its time variation dA/dt is also smaller. Further, it 
is known that the deviation of temperature in the observed cosmic microwave back- 
ground (CMB) radiation is of the order of 10^^; hence, we may use the following, 

« 10"^ (37) 

( LiA\ 

\A dt ) 

It is, therefore, currently difficult to detect the cosmological contributions attributed 
not only to the RW model but also to the inhomogeneity of the universe. The inhomo- 
geneity in the background cosmological matter distribution does not have a detectable 
effect and hence cannot explain the observed dAU/dt. 



5 Conclusions 

We investigated the cosmological influence according to the LTB model on a grav- 
itationally bound local system such as the solar system. We focused on planetary 
motion and obtained the leading-order dynamical perturbation generated from the 
LTB spacetime. The obtained dynamical perturbations, especially ( l34l l. contains the 
contribution attributed to the RW model and also the correction terms attributed to 
the LTB model. Moreover, we applied the obtained results to the secular increase in 
the astronomical unit 1491 and confirmed that the effect of the inhomogeneity of the 
universe does not provide an explanation for dAU/dt. 
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In spite of several attempts (see section[T]i, the origin of dAU / dt is far from clear. 
It is now pointed out that the most plausible reason of dAV/dt is due to either the 
lack of calibrations in the internal delays of radio signals within spacecrafts or the 
complication of the modeling of solar corona. However, none of the explanations have 
thus far been successful; this issue should hence be explored using all possibilities. A 
re-analysis of dAU/dt incorporating new data sets is also expected. 

Since the astronomical unit, as expressed in ([2), is currently determined from the 
arrival time measurement of radar signals, we must investigate this problem in terms 
of light/signal propagation. To this end, we need to construct a cosmological model 
that combines the LTB model with a gravitating body. Of course, it may be difficult to 
detect the cosmological effect in the solar system. Nonetheless, since the theoretical 
discussions pertaining to this issue are still unresolved, it is, from theoretical point of 
view, interesting to develop a rigorous physical model that matches the gravitational 
bound local system and the cosmological models and clarify several assertions. 

Acknowledgements We would like to thank the referees for comments and suggestions. We also ac- 
knowledge Prof. G. A. Krasinsky for providing information and comments regarding the AU issue. This 
work was partially supported by the Ministry of Education, Science, Sports and Culture, Grant-in-Aid, No. 
21740193. 

References 

1. Standish, E. M., JPL Interoffice Memorandum, 312N, 03 (2003) 

2. Pitjeva, E. V., Solar Syst. Res., 39, 176 (2005) 

3. Bretagnon, P, Francou, G., A&A, 202, 309 (1988) 

4. Fienga, A., Manche, M., Laskar, J., Gastineau, M., A&A, 477, 315 (2008) 

5. Soffel, M. et al, AJ, 126, 2687 (2003) 

6. Will, C. M., Theory and Experiment in Gravitational Physics, Cambridge Univ. Press (1993) 

7. Will, C. M., Living Reviews in Relativity, |http://relativity.livingreviews.org/Articles/lrr-2006-3/| (2006) 

8. Danzmann, K., Adv. Space Res., 25, 1129 (2000) 

9. Turyshev, S. G., Shao, M., Nordtvedt, K., Class. Quant. Grav., 21, 2773 (2004) 

10. Ni, W.T., Int. J. Mod. Phys. D, 17, 921 (2008) 

11. Damour, T, Soffel, M., Xu, C, Phys. Rev. D, 43, 3273 (1991) 

12. Damour, T, Soffel, M., Xu, C, Phys. Rev. D, 45, 1017 (1992) 

13. Damour, T, Soffel, M., Xu, C, Phys. Rev. D, 47, 3124 (1993) 

14. Damour, T, Soffel, M., Xu, C, Phys. Rev. D, 49, 618 (1994) 

15. Bramberg, V., Kopeikin, S., Nuovo Cimento B, 103, 63 (1989) 

16. Brumberg, V., Kopeikin, S., Reference frames in astronomy and geophysics, 1 15 (1989) 

17. Perlmutter, S. etal, ApJ, 517, 565 (1999) 



Title Suppressed Due to Excessive Length 



11 



18. McVittie, G. C, MNRAS, 93, 325 (1933) 

19. Jiimefelt, G., Ann. Acad. Soc. Sci. Fennicae, A45, 3 (1940) 

20. Jamefelt, G., Arkiv Matem. Astron. Fys., 27, 1 (1940) 

21. Jamefelt, G., Ann. Acad. Soc. Sci. Fennicae, A45, 12 (1942) 

22. Einstein, A., Straus, E. G., Rev. Mod. Phys., 17, 120 (1945) 

23. Einstein, A., Straus, E. G., Rev. Mod. Phys., 18, 148 (1946) 

24. Schilcking, E., Zeitschrift fiir Physik, 137, 595 (1954) 

25. Noerdlinger, P D., Petrosian, V., ApJ, 168, 1 (1972) 

26. Cooperstock, EI., Earaoni, V., VoUick, D. N., ApJ, 503, 61 (1998) 

27. Khoner, S. A., Soffel, M. H., Proc. of the Symposium "The Three-Dimensional Universe with Gaia", 
305 (2004) 

28. Earaoni, V., Jacques, A., Phys. Rev. D, 76, id 063510 (2006) 

29. Sereno, M., Jetzer, P, Phys. Rev. D, 75, id 064031 (2007) 

30. Adkins, G. S., McDonnell, J., Phys. Rev. D, 75, id 082001 (2007) 

31. Kopeikin, S. M., AIP Conference Proc. "New Trends in Astrodynamics and Applications IIP', 886, 
268 (2007) 

32. Carrera, M., Giulini, D., Rev. Mod. Phys., 82, 169 (2010) 

33. Tomita, K., ApJ, 529, 26 (2000a) 

34. Tomita, K., ApJ, 529, 38 (2000b) 

35. Tomita, K., MNRAS, 326, 287 (2001a) 

36. Tomita, K., Prog. Theo. Phys., 105, 419 (2001b) 

37. Tomita, K., Prog. Theo. Phys., 106, 929 (2001c) 

38. Iguchi, H., Nakamura, T., Nakao, K., Prog. Theor. Phys., 108, 809 (2002) 

39. Lematoe, G., Ann. Soc. Sci. Bruxelles A, 53, 51 (1933) 

40. Tohnan, R. C, Proc. Nat. Acad. Sci., 20, 410 (1934) 

41. Bondi, H., MNRAS, 107, 410 (1947) 

42. Plebahski, J., Krasihski, A., An Introduction to General Relativity and Cosmology, Cambridge Univ. 
Press (2006) 

43. Kasai, M., Prog. Theor. Phys., 117, 1067 (2007) 

44. Gautreau, R., Phys. Rev. D, 29, 198 (1984) 

45. Mashhoon, B., Mobed, N., Singh, D., Class. Quant. Grav., 24, 5031 (2007) 

46. Krasiliski, A., Inhomogeneous Cosmological Models, Cambridge Univ. Press (1997) 

47. Anderson, J. D., Laing, P. A., Lau, E. L., Nieto, M. M., Turyshev, S. G., Phys. Rev. Let., 81, 2858 
(1998) 

48. Anderson, J. D., Campbell, J. K., Ekelund, J. E., ElUs, J., Jordan, J. F., Phys. Rev. Lett., 100, id 091102 

(2008) 

49. Krasinsky, G. A., Brumberg, V. A., Celest. Mech. Dyn. Astm., 90, 267 (2004) 

50. lorio, L., AJ, 137, 3715 (2009) 

51. Lammerzahl, C, Preuss, O., Dittus, H., Lasers, Clocks and Drag-Free Control: Exploration of Rela- 
tivistic Gravity in Space eds H. Dittus, C. Lammerzahl, and S. G. Turyshev, Springer, Berlin, 75 (2008) 

52. Standish, E. M., Proc. lAU CoUoq., 196, 163 (2005) 

53. Pitjeva, E. V., Private Communication (2009) 

54. Krasinsky, G. A., Private Communication (2007) 



12 



Hideyoshi Arakida 



55. Arakida, H., New Astron., 14, 264 (2009) 

56. Noerdlinger, P. D., arXiv:0801.3807 (2008) 

57. Arakida, H., Adv. Space Res., 45, 1007 (2010) 

58. Weinberg, S., Gravitation and Cosmology, John Wiley & Sons, New York, Chichester, Brisbane, 
Toronto, Singapore (1972) 

59. Soffel, M. H., Relativity in Astrometry, Celestial Mechanics and Geodesy, Springer Verlag, Berlin, 
Heidelberg, New York, London, Paris, Tokyo (1989) 

60. Robertson, H. R, Phil. Mag., 7, 845 (1928) 

61. Irvine, W. M., Ann. Phys. (N.Y), 32, 322 (1965) 

62. Kodama, H., Sasaki, M., Prog. Theo. Phys. Suppl., 78, 1 (1984) 

63. Tomita, K., Prog. Theo. Phys., 85, 1041 (1991) 

64. Mukhanov, V. R, Feldman, H. A., Brandenberger, R. H., Phys. Rep., 215, 203 (1992) 

65. Shibata, M., Asada, H., Prog. Theo. Phys., 94, 1 1 (1995) 

66. Dodelson, S., Modem Cosmology, Academic Press (2003) 

67. Partovi, M. H., Mashhoon, B., ApJ, 276, 2 (1984) 



